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Abstract 

Properties of Jost and dual Jost solutions of the heat equation, k) and ^{x, k), 
in the case of a pure solitonic potential are studied in detail. We describe their analytical 
properties on the spectral parameter k and their asymptotic behavior on the a;-plane 
and we show that the values of e~''^^{x,'k.) and the residua of e''^*I'(2;, k) at special 
discrete values of k are bounded functions of a: in a polygonal region of the g-plane. 
Correspondingly, we deduce that the extended version L{q) of the heat operator with a 
pure solitonic potential has left and right annihilators for q belonging to these polygonal 
regions. 

00 ■ 

\o 

^ ' 1 Introduction 
in 

' The Kadomtsev-Petviashvili equation in its version called KPII 

{Ut - 6uUx^ + Uxixixi)xi = -3Ux2X2, (1-1) 

. , , where u = u{x,t), x = (xi,X2) and subscripts xi, X2 and t denote partial derivatives, is 

r> ■ a (2+l)-diniensional generalization of the celebrated Korteweg-de Vries (KdV) equation. 

I There are two inequivalent versions of the KP equations, corresponding to the two choices 

± for the sign in the rhs of which are referred to, respectively, as KPI and KPII. The 

KP equations, originally derived as a model for small-amplitude, long-wavelength, weakly 
two-dimensional waves in a weakly dispersive medium [T] , were known to be integrable since 
the beginning of the 1970s [31 [5], and can be considered as prototypical (2+l)-dimensional 
integrable equations. 

The KPII equation is integrable via its association to the operator 

Cix,dx)^-dx,+dl~u{x), (1.2) 

which defines the well known equation of heat conduction, or heat equation for short, since 
it can be expressed as compatibility condition [£,7^ =0 of the Lax pair C and T , where 
7" is given by 

T{x, dx,dt) ^dt + 49^j - 6udxi ~ iux^ - Sd'^Ux^. (1.3) 

The spectral theory of the operator (|1.2p was developed in [H [5J [HI [7] in the case of 
a real potential u{x) rapidly decaying at spatial infinity, which, however, is not the most 
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interesting case, since the KPII equation was just proposed in [T] in order to deal with two 
dimensional weak transverse perturbation of the one soliton solution of the KdV. In fact, if 
Uj^(t,Xi) obeys KdV, then u{t, Xi,X2) = Ui{t,Xi + /xxj ~ solves KPII for an arbitrary 

constant /i G R. In particular, KPII admits a one soliton solution of the form 



u{x, t) 



(k1 - K2f 



■sech 



-xi 



'-X2 - 2{kI - Kl)t 



(1.4) 



where ni and K2 are real, arbitrary constants. 

A spectral theory of the heat operator ()1.2p that also includes solitons has to be built. 
In the case of a potential u{x) = uo{x) rapidly decaying at spatial infinity, according to [H 
O O [7], the main tools in building the spectral theory of the operator (|3.ip . as in the 
one dimensional case, are the integral equations, whose solutions define the related Jost 
solutions. However, if the potential u{x) does not decay at spatial infinity, as it is the case 
when line soliton solutions are considered, the integral equations of the decaying case are 
ill-defined and one needs a more general approach. In solving the analogous problem for the 
nonstationary Schrodinger operator, associated to the KPI equation, the extended resolvent 
approach was introduced [8] . Accordingly, a spectral theory of the KPII equation that also 
includes solitons has to be investigated using the resolvent approach. In this framework it 
was possible to develop the inverse scattering transform for a solution describing one soliton 
on a generic background ^ , and to study the existence of the (extended) resolvent for (some) 
multisoliton solutions [TU]. 

However, the general case of 7V-solitons is still open. Following [8], the first step in 
building the inverse scattering for this case lies in building a Green's function G{x, x' , k) of 
the heat operator (|1.2|) corresponding to the pure A^-soliton potential u{x) such that 



.jH^'i-^i)+^'i^'2-^2)G{x,x',k) 



(1.5) 



is bounded for any x and any k G C. We expect this Green's function to be bilinear in 
the Jost solutions <l'(x,k) and ^(x, k) of the heat operator and its dual and, therefore, we 
need, first, to study the general properties of these solutions. Once performed this study 
and obtained the Green's function G, the Jost solution $ for a potential given as a sum of 
two terms, 

u{x) — u(x) + v! (x\ (1-6) 

where u{x) is a pure A^-soliton potential and it'(x) an arbitrary bidimensional decaying 
smooth perturbation, can be derived as solution of the integral equation 



$ = $ + Gu'*, 



(1.7) 



that is well defined since, thanks to u', it has a good behaving kernel Gv! . Analogously, for 
the Jost solution ^I^ of the dual heat operator. Then, having at hands these building blocks, 
$ and one can assemble the inverse scattering theory. 

The paper in organized as follows. In Sec. 2 we consider the multisoliton potential 
showing at large spaces iV;, "ingoing" rays and Na "outgoing" rays obtained in its general 
form in |llj . and we reformulate it in terms of r- functions as in the review paper |12j and 
in ^13]. In Sec. 3 we consider the Jost solutions of the heat operator and their dual, already 
obtained in [llj , but they also expressed in terms of r- functions as in [13^ . In view of getting 
the asymptotic behavior of the Jost solutions, in Sec. 4 the asymptotic behavior at large x of 
the multisoliton potential, studied in details in jl3] , is reviewed showing that the round angle 
at the origin can be divided in A/" = iV^ + iVb angular sectors such that on their bordering 
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rays the potential has constant soliton like behavior, while along directions inside the sectors 
the potential has an exponential decaying behavior, which is explicitly given. In Sec. 5, by 
using the results obtained for the potential, also the asymptotic behavior at large x of the 
Jost and dual Jost solutions and their values or, accordingly, residues at the discrete points 
of the spectrum are obtained. In Sec. 6 by using the results obtained in the previous section 
we show that the values of e~'^^^{x,'k) and the residue of e''^^'(x,k) at special discrete 
values of k are bounded functions of a; in a polygonal region of the g-plane included inside 
the parabola q2 = qf- Correspondingly, we conclude that the extended version L{q) of the 
heat operator with a pure solitonic potential has left and right annihilators for q belonging 
to these polygonal regions. It follows that the extended resolvent M{q), inverse of L{q), 
does not exists for q belonging to these polygons. However, the Green's function G{x, x', k), 
needed, as indicated above, for building the inverse scattering theory in the case of KPII 
solutions with solitons, can be obtained by a reduction procedure involving q values outside 
the parabola. 

2 Multisoliton potentials and their properties 
2.1 Main notations 

Soliton potentials are labeled by the two numbers (topological charges), Na and Nb, that 
obey condition 

Na,Nb>l, (2.1) 

Let 

N'^Na+ Nb, (2.2) 
so that N >2. We introduce the TV real parameters 

Kl < K2 < . . . < nu, (2.3) 

and the functions 

K^{x) ^ K,^x-^+ K^x^, n = l,...,A/'. (2.4) 

Let 

e^(-)=diag{e^"(-)}^^i, (2.5) 

be a diagonal J\f x TV- matrix, T) &J\f xNb constant matrix with at least two nonzero maximal 
minors, and let V be an "incomplete Vandermonde matrix," i.e., the Nb x A/'-matrix 



V 



/I ... 1 



(2.6) 



Then, the soliton potential is given by 

u{x) ^ -2dl^logTix), (2.7) 
where the T-function can be expressed as 

t(x) =det(Ve^(^)2?). (2.8) 

See the review paper |12| . references therein, and |13) where the same notations have been 
used. 
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There exists (see [Tl], [TU], and [T3]) a dual representation for the potential in terms of 
the T-function 

T'(a;) =det(l?' 6-^(^)7 V') , (2.9) 

where I?' is a constant Na x A/'-matrix that like the matrix D has at least two nonzero 
maximal minors and that is orthogonal to the matrix P in the sense that 

P'P^O, (2.10) 

being the zero in the r.h.s. a Na x A^b-nratrix, and where V' is the TV x A'^^- matrix 

/I ... 

: : I , (2.11) 

V 1 ... 

and 7 the constant, diagonal, real Af x A/'-matrix 



AT 

7 = diag{7„}^^l, Jn ^ II i'^n - Kn')-\ (2.12) 

n' = l 
n' 

In order to study the properties of the potential and the Jost solutions, it is convenient 
to use an explicit representation for the determinants. By using the Binet-Cauchy formula 
for the determinant of a product of matrices we get 

l<ni<n2<---<nNj<7V 1=1 

with 

= F(K„^,...,K„^JX»(ni,...,nArJ, (2.14) 

where we used notation 

V(Kl,...,K^f) = Y\_ (l^n-Km), (2-15) 
l<m<n<Af 

for the Vandermonde determinant and notation 

: : , (2.16) 

for the maximal minors of the matrix 2?. Notice that the coefficients fni,n2....,nN^ ^-re invari- 
ant under permutations of the indices. 

In what follows it is convenient to consider indices n, rii, . . . , rijVi, of parameters k„, 
coefficients fni,n2,...,nN^j 6tc, to be running modTV throughout Z, i.e., let 

n n{modJ\f). (2.17) 

From (|2.13p it follows directly that condition 

/«!,...,«„, > 0, for aU 1 < ni < na < • • ■ < njv^ < A/", (2.18) 
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is sufficient (see [TS]) for the regularity of the potential u{x), i.e., for the absence of zeros of 
t{x) on the x-plane. Thanks to (|2.3I) . (|2.15p and (|2.14p this condition is equivalent to the 
condition that all maximal minors of the matrix D are non negative. In |13| it was mentioned 
that condition (|2.18p is also necessary for the regularity of a potential under evolution with 
respect to an arbitrary number of higher times of the KP hierarchy. In |15| it is suggested 
to decompose soliton solutions of KPII into subclasses, associated to the Schubert cells on 
Grassmanian. It is proved there that condition (|2.18p is necessary for the regularity of the 
all solutions associated to a cell. However, the problem of finding necessary conditions for 
the regularity of the many soliton solution under evolution with respect to KPII only is still 
open. 

Under condition (|2.17l) and thanks to (|2.3|) we have the following lemma (see [TB]) that 
we use below. 

Lemma 2.1 Let 

gi,m,n = {kI - Km){Kn - '«n+ArJ(K; + Km - Kn - Kn+AfJ- (2-19) 

Then 

9i,m,n^^ for any n G Z, I = n, . . . ,n -\- Njj, m ^ n + Nj,, . . . ,J\f +n, (2.20) 

and equality takes place only when I and m independently take values n or n + Nb by modJ\f. 

Matrices V and T>' obey rather interesting properties, see [131 [in]- In particular, products 
V and V' 'D'\ where f denotes Hermitian conjugation of matrices (in fact, transposition 
here), are positive and then invertible, so there exist matrices (see |17j ) 

(2?)'"^^ = (2?t2?)-i2?t^ (py^^^ ^2?'t(2?'2?'1')-i, (2.21) 

that are, respectively, the left inverse of the matrix T> and the right inverse of the matrix 
V\ i.e., 

{vY-'^V^En„ V'{vf-'^ ^En^, (2.22) 

where E^^ and Ejy^ are, respectively, the Na x Na and N), x Ni, identity matrices. Products 
of these matrices in the opposite order give the real self-adjoint M x TV-matrices 

P = V{V)'-^^^ =V{V^V)-^{V)\ (2.23) 

P' = {V')^'^^V' = {V')\V'V'^)-^V', (2.24) 

which are orthogonal projectors, i.e., 

P^ = P^ {P'f^P\ pp'^Q^p'p^ (2.25) 

and complementary in the sense that 

P + P'^Eu- (2.26) 

Matrices D and T)' are not in one-to-one correspondence with potential u{x). Indeed, by 
(|2.7p . (|2.8p . ()2.9|) one gets that the potential is invariant under substitutions 

V^Vv, V'^v'V', (2.27) 

where v and v' are arbitrary constant, nonsingular Nh x N^,- and Na x iVo-matrices, corre- 
spondingly. Thus under condition (|2.3I) (A^a, A^b)-soliton potential is parameterized by the 
point of the Grassmanian GTN^,j\f, if representation (j2.8|) is used, or by the point of the 
Grassmanian Grjv„Xi if representation p.9p is used. 
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Figure 1: Sector (T„ corresponds to Na + 2 < n < JV — 1. 

2.2 Asymptotic behavior of the r-function and potential u{x) 

Here we shortly report (see [IB] for details) the main asymptotic properties of the function 
t{x) and potential u{x) at large space, which are necessary to know in order to derive the 
asymptotic properties of the Jost and dual Jost solutions and of their discrete values. 

The asymptotic of t{x) is determined by the interrelations between functions Kn{x) for 
different n, that is, due to (|2.4p . by the differences 

Ki{x) - Km{x) = [ki - Km){xi + (k/ + Km)a;2), for any l,m eZ, (2.28) 

which are linear with respect to the space variables. Then, the asymptotic behavior must 
have sectorial structure on the x-plane and in order to describe this structure we divide the 
a;-plane in sectors. 

We introduce on the a;-plane rays r„ given by 

rn = {x : Xi + {Kn+N^ + Kn)x2 = 0, {Kn+N^ ~ Kn)x2 < 0), 71=1,..., A/". (2.29) 

They intersect the lines X2 = ±1, respectively, at the points 

('tn+ATfc + Kn, -1) for n = 1, . . . ,7Va 

{-Kn+N,- Kn,l) ioi Tl ^ Na + I , . . . , M ■ (2.30) 

Therefore, increasing n from n = 1 to n = Na the ray rotates anticlockwise in the lower half 
x-plane, crosses the positive part of the a;i-axis in coming to n = Na + 1 and, then, rotates 
anticlockwise in the upper half a;-plane up to n = A/" and, finally, crossing the negative part 
of the xi-axis, for n = Af + 1 comes back to the ray ri{x) thanks to (|2.17p . see Fig.[T] 

Let us assume that some rays, say, and r„, where for definiteness m < n, are parallel. 
By (|2.29p this means that K^n+Nt + i^m = Kn+Nt + >^n, i.e., that k„ - Km = ^m+Ni - ^n+N^, 
where the l.h.s. is positive thanks to (|2.3I) . Then, because of (|2.17p it is easy to see that 
the r.h.s. can be positive only if 1 < m < TVa < n < A/", so that by (|2.30l) ray is in 
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the bottom halfplane and r„ in the upper one. We see that for generic values of the k„'s 
and when Na ^ Nh the rays in the upper and bottom half planes cannot be parallel, while 
this is possible for a special choice of k„'s. On the contrary, in the case Na = Nb all pairs 
r„ and r„-|_Ar^, n = 1, . . . ,Na, and only these pairs give parallel rays. In the special case 
Na = Nh = 1 we get two rays producing the straight line xi + {ki + K2)x2 = that divides 
the cc-plane in two halfplanes. 

In the same way we introduce sectors (t„ , which are subsets of the x- plane characterized 

as 

(7n = {x : Kn-i{x) < i^„+Ar^_i(x) and Kn{x) > Kn+NA^)}^ for n = 1, ... ,7V. (2.31) 

Since (I2.28P and the discussion above it follows that the ct„ 's are sharp (for J\f > 2) angular 
sectors with vertices at the origin of the coordinates bounded from the right (looking from 
the origin) by the ray Vn-i and from the left (looking from the origin) by the ray r„ and 
that the sectors (T„ with the increasing of n are ordered anticlockwise, starting "from the 
left" with the sector ci, that includes the negative part of the xi-axis, and, then, with the 
sectors an {n — 2, . . . , Na) in the bottom half-plane, the sector crjv„+i "to the right," that 
includes the positive part of the xi-axis, and the sectors tT„ (n = Na + 2, . . . ,Af) on the 
upper half-plane, finishing with the sector aj\f tangent to the sector cti, covering in this way 
the whole x-plane with the exception of the bordering rays r„, see Fig.[TJ 

Therefore, the sectors cr„ define a A/'-fold discretization of the round angle at the origin, 
with the integer n (mod J\f) playing the role of a discrete angular variable. It is clear that in 
the study of the asymptotic behavior of the r-function we consider the x-plane as a vector 
space, since the finite part of x when a; ^ cx) is irrelevant. 

For determining the directions of rays r„ and sectors cr„, we introduce the vectors 

2/n = (S^Aff, - «n - Kn+ATt), n = 1 , . . . , TV, (2.32) 

and use them for giving the following definition. 

Definition 2.1 We say that x oo along the ray Vn if x ^ oo and there exists such 
a — > -l-cxD that X — ayn is bounded. This will be denoted as x — ^ oo. 

We say that x ^ oo in the sector cr„, n = 1, . . . ,Af, if x ^ oo and there exist such 
a — > +00 and /3 — > +00 that x — ayn-i — f3yn is bounded. This will be denoted as x — ^ 00. 

Notice that for x 00 

Knix) - Kn+N^{x) IS bouuded and (i^n+Nb - i^n)x2 -00, (2.33) 

and for x 00 

Kn+Nt-l{x) - Kn-l{x) +CX1, Kn{x) - K^+nA^) ^ +°°^ (2.34) 
as follows directly from the definition. In fact we have a more general statement. 
Lemma 2.2 Let Na, Ni, > 1, and n E Z arbitrary. Then 

1. if X oo we have that Ki{x) — Km{x) — > +oo or bounded for any I ~ n, . . . ,n + Nb 
andm = n+N^, . . . ,Af+n, where boundedness takes place if and only if(l,m) = (71,71), 
(71, 71 + Nb), (n -I- Nb, n), {n + Nb,n + Nb); 

2. if X 00 we have that Ki{x) — K„i{x) — > +cxd for any I = n, . . . ,n + Nb — 1 and 
m = n + Nb,...,Af+n- 1; 
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where summation of indices is always understood mod Af. 



Thanks to the above Lemma [2.21 the asymptotics of the r- function is given by the fol- 
lowing Theorem (see [H]). 

Theorem 2.3 // condition !i2.18]} is satisfied, the asymptotic expansion of t{x) for x — )■ oo 
for any n is given by 



^{x) = (z„ + z„+ie^"^+"(-)-^"(-) + 0(1)) exp V K,{x) , (2.35) 




X — > oo : 

x^oo: t(:e) = (z„ + o(l))exp f ^ if/(x)j , (2.36) 
where notation 

Zn = /«,«+!,. ..,n+W,-l = V{k^, K.^_^^^_^)V{n, . . . , U + Nj, - I) , (2.37) 

was introduced for any n = 1, . . . ,7V. 

As we see from (|2.35p and (|2.36p the leading asymptotic term along the ray direction 
r„ is given by the sum of the leading terms obtained for x oo and x oo. Since 
the exponential factor cancels out when this expansion (I2.35|) of t{x) is inserted in (12.71) . 
the factor in parenthesis gives the ray behavior of the potential u(x) at infinity. Explicitly, 
taking into account (|2.17p we get Na asymptotic rays in the bottom half plane: X2 — oo, 
xi + {ktl + Kn+Nb)x2 boundcd. Along these rays the potential behaves as 

M(x) = -292^1og(z„ + z„+ie^"+"<.(^)-^"(^)), n = l,...,Na. (2.38) 

In the same way, the potential u{x) has A^^, asymptotic rays in the upper half-plane: X2 
-|-(X), Xi -\- (Kji 4" /^n+ Na )x2 bounded. Along these rays it behaves as 

u{x) - -2dl^ log(z„+jv„ + ^n+iv^+ie^"^")-^"^"'''"'), n^l,...,m. (2.39) 

Asymptotic behavior inside sectors (t„ is given by (j2.36l) where the only a;-dependent term 
is the exponential factor. Taking into account its linear dependence on x and (I2.18P we get 
that u{x) decays (exponentially) in all directions inside all sectors, i.e., on the whole x-plane 
with exception to the rays ri, . . . , r^f. 

Thanks to condition p.l7p coefficients z„ defined in ()2.37p are a special subset of coeffi- 
cients /ni,....n„j_ in (|2.14p . Theorem 12.31 shows that the behavior of the potential at large x 
is determined by the coefficients z„ only. Therefore, condition 

z„>0, n = l,...,7V, (2.40) 

is a sufficient condition for the regularity of the potential at large x. If the stronger condition 
(|2.18p is not satisfied, singularities can appear, but only in a finite region of the x-plane. 
Taking into account the special role played by coefficients z„ let us introduce matrices 



n^l,...,M, (2.41) 
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so that by ([^1^ det w„ = V{n, . . . ,n + Nb - 1) and by (P37| 



Zn = V"(k„,...^k„+jv^_i ) det Vn, (2.42) 
that is different from zero thanks to (I2.40p . Then, we can perform the special permutation 

7T„{l)^Na + l + l-n, l^l,...,JV, (2.43) 
that shifts matrix w„ in the bottom part of matrix J). Thus we get the representations 

V = n,,(^ ) ^' = -^mH^ (2-44) 

where the Na x iV^-niatrix is defined in analogy. Here tt„ denotes matrix of permutation 
(|2.43l) . TT^j its Hermitially conjugate, E]\[^ and E^^ are, correspondingly, Na x A'^a and Nt x A^f, 
unity matrices. The Na x A'^^ rectangular submatrix d(„) is defined by any of the equalities 
in (|2.44l) and it is obvious that, thanks to their block structure, matrices V and V' obey 
(|2.10p . Representation (|2.44p is convenient for studying the asymptotic properties of the 
Jost solutions. 



3 Jost solutions 

Since the heat operator is not self-dual, one has to consider simultaneously its dual C^ix, dx) = 
+ — u{x) and, then, introduce the Jost solution $(a;, k) and the dual Jost solution 
\['(a;,k) obeying equations 

£(a;,9,)$(a;,k) =0, /:''(x, a,)*(a;, k) = 0, (3.1) 

where k is an arbitrary complex variable, playing the role of a spectral parameter. The 
reality of the potential u{x), that we always assume here, is equivalent to the conjugation 
properties 

$(a;,k) = $(a;, -k), ^'(a;,k) = ^-(2;, -k). (3.2) 
The T-function representations for the Jost solutions were derived in |13j and are given as 

$(x,k) ^ "'^(^' '^) e-'kxi-k^..^ ^(.T,k) = ^'i^(^' '^) e'k^i+k^^2^ (3 3) 

t{x) t{x) 

with 

rcE.(x,k) =dct(Ve-^(^)(K + ik)P), T^ix,k) = det(Ve-^(^)(K + z k)-^ D), (3.4) 
where k + i k denotes the diagonal M x A/'-matrix 

/t + i k = diag{K]^ + i k, . . . , + i k} (3.5) 

and analogously for the matrix (k + ik)^^. By this definitions e^^^^^^ ^2<i>(a::,k) is a 
polynomial with respect to k of order Ni, and g-'^'^^i-'^ ^^'^/{x, k) is a meromorphic function 
of k that becomes a polynomial of order Na after multiplication by Jl^i i'^n + « k). In other 
words, $(a;, k) is an entire function of k and "^{x, k) a meromorphic function with poles at 
points k = iKn, n = 1, . . . ,7V. Introducing the discrete values of $(a;, k) at these points as 
a A/"- row 

$(a;, in) = {$(x, in^), ^{x, iK^f)}, (3.6) 
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and the residuals of ^{x, k) at these points 

^-^Jx) = res ^'(a;,k), (3.7) 

as a A/"- column 

^^ix)^{^^,ix),...,^^^ix)}^, (3.8) 
we have the fohowing relations (see [T3]) 

^{x,iK)V = 0, 2?'*„(a;)=0. (3.9) 

Thanks to (|2.23p and (|2.24p relations p.9p can be written equivalently in the form ^{x, iK)P = 
0, P'-^i^{x) = 0, so that by (P?^ and the last equality in (P?^ 

^{x, in) = $(a;, iK){P + P') ^ $(a;, *^(a;) = (P + P')*K(a;) = -P*«:(a;), (3.10) 

and, then, by (|2.25l) for any x and x' 

$(a;,iK)^'«(x') = ^$(a;,iKn)*K„(a:') =0. (3.11) 

n=l 

This means that the product $(a;, k)vE'(a;', k) of the Jost solutions obeys the well known 
Hirota bilinear identity [18] for the Beiker-Akhiezer solutions, if the contour of integration 
surrounds all points k = iK„. 
Let us introduce 

ip{x) ^ <^>{x,iK)V'^~'^\ ^{x) =V^-^^ "ii^ix), (3.12) 
where notation (|2.2ip was used. By this definition, we have 

(p{x) ^ {(pi{x),...,(pN^{x)), Na-TOW, (3.13) 

V'(x) = (?/^i(a;), . . . ,?/'jVi,(a;))'^, iVfe-column, (3.14) 
and, by (|2.23l) and (I2.24p . we get representations 

<i>{x,iK) = (p{x)V', ■$^ = Vip{x). (3.15) 

Thus we constructed Na solutions and Nt dual solutions that parameterize discrete values 
of the Jost and dual Jost solutions (that have TV components). It is clear that (f{x) and 
ipix) are not invariant with respect to the redefinition of matrices V and D mentioned in 
(j2.27p . while, in contrast, combinations T>v^^ and v'^ D' are invariant with respect to 
transformation (|2.27l) . as follows from (|2.44p . 

It is necessary to mention that with the above definitions p.3p the Jost solutions have 
at large k the asymptotics 

lim (ik)"^''e**'^i+'''^^$(x,k) = 1, lim (i k)^''e-*''^i-''' ^^^'(x, k) = 1, (3.16) 

k— >oo k— )-oo 

and the potential is reconstructed as 

u{x) = -2 lim {ik)-^''+^d^Ae'^^'+^' ^^^{x,k)) = 

k^oo 

= 2 lim {ik)^''+^d^Ae-'^'=''^''-'^^{x,k)). (3.17) 

k^oo 
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4 Asymptotics of the Jost solutions 



4.1 Asymptotics of the Jost solutions $(a;, k) and k) for a generic 

keC 

Since, as was already noted in [13 , the functions T$(x,k) and T^(x,k) defined in (j3.4p can 
be obtained from the function t{x) by means, respectively, of the special Miwa shifts [18] 
^K{x) _^ e^^^)(K + ik) and e^*^^-* e^^^)(K + i k)~^, their asymptotic behavior follows 
trivially from Theorem 12.31 for t{x) if k 7^ iHn for all n. Therefore, from equation (I3.3P we 
get that the Jost solutions for x 00 inside an arbitrary angular sector an, n = 1, . . . , Af, 
have the following asymptotic behaviors: 

$(x,k) =6-^''"^-'''"^ n (Kj+Zk) + ..., (4.1) 

j=n 
n+Ni-1 

^{x,-k) =e'^^'+^'^' Yl (k, +zk)-i + ..., (4.2) 

j=n 

were summation of indices is by mod TV and were dots denote weaker terms. Asymptotic 
behavior along the ray, x 00 also easily follows from Theorem 12.31 and we omit it here, 
since in what follows we are interested only in the exponential behavior of the leading terms, 
while in the ray directions we have two leading terms with the same asymptotic behavior 
(see (|2.35p ). To emphasize this we introduce the "closed" sectors a^- In other words we 
introduce: 

Definition 4.1 We say that x 00 along the direction of the closed sector a^, n = 
1, . . . ,J\f, if X 00 and there exist non negative a and l3 such that x — ayn-i — Pyn is 

bounded when a + /3 ~> 00. This will be denoted as x — ^ 00. 

According to this definition, asymptotics (|4.ip and (|4.2I) are valid for x 00. 

In what follows it is enough to study the asymptotic behavior of only one of the Jost 
solutions, $(a;,k) or \E'(a;,k). Indeed, thanks to p.3p . the dual representations for r in (|2.8|) 
and (|2.9|) and the analogous ones for r$ and (see Remark 5.1 in [l3]), one can derive that 
the Jost solutions ^{x,k) and ^'(a:, k) are related by the equation 

AT 

*(x,k) = [$(-a;,k)] [](K„ + zk)-\ (4.3) 

n=l 

where square brackets in the r.h.s. means that one has to perform in $(— x, k), as defined 
in p.4p . the substitutions 

Na^Nb, V^jV'^, (4.4) 

where matrices V ' and 7 are defined in (|2.10p and (|2.12p and superscript T denotes matrix 
transposition. It is easy to see that asymptotics (|4.ip and (|4.2p obey this property, if one 
notice that, thanks to Definitions 12. H substitution (|4.4p leads to the following relations for 
the sectors: 

[fn] = (Jn+Na^ Wn+Na] = m (4.5) 

SO that asymptotic behavior of '^{x,\<i) in sector cr„ is given by the asymptotic behavior of 
$(— a;,k) in sector (7„+Ar^, once performed the transformations (|4.4p . 



11 



When k = iKm, for some m, some terms in r$(2:,k) (see p.4p ) are absent. In the case 
when such terms are coefScients of the leading exponents we get zero asymptotics and exact 
behavior of the Jost solution cannot be derived from (14. ip . The same is valid for residuals 
of ^(x, k). Thus we have to consider asymptotic behavior at these values of k separately. 



4.2 Asymptotics of the discrete values of the Jost solutions 

In order to find the asymptotic behavior of the discrete values of the Jost solutions, we 
notice that according to p.3p . p.4p and definition (|2.4p . they equal 

/ • \ res T^(a;,k) 

= !M^L!^e^™W, *,^^(^) ^ ^^^Hl^ ^-K^(^)^ (4.6) 

t[x) t[x) 

wherem = 1,...,M. We consider the asymptotic of $(a;, ik™), since the behavior of ^P^^ (a;) 
can be obtained from (14.31) that for these discrete values means 



(x) = ^(-l)^7rn [H-X, ll^rn)] , (4.7) 

where we used notation (|2.12p and where the square brackets denotes operation (14. 4p . 
Let us fix some n = 1, . . . , A/" and let x — ^ oo. From (j4.1l) we have 

am\x), m = n, . . . ,n + Nb — I, 



I Pm e TO — n + iVb,...,ri+7V— i, 

where constants pm"* equal 



n («. -«™). (4-9) 

and where functions a™ (x) , that give the leading asymptotic behavior for the corresponding 
values of m, should be determined (thanks to (|4.1I) we know only that they tend to zero 

when X — ^ oo). 

Using the permutation introduced in (j2.43l) we write 

$(x,iK)7r„ = [^{x,iKm+n+Nt-l))^,^i = 

Now by ^^M^ we get from ([53)1 

where we used that matrices Vn are invertible. Inserting these relation in (|4.10l) and again 
using (I2.44P we get 

$(x,z«) = (pi'^^^e^-". (-),..., pi"^^_,e^"+--(-)).:r^ 2?' + ..., (4.11) 
*.(x) = -zP«„-i((p("))-^e-^"(^), . . . , (pit^,_0~'e-^— ^(^))% . . . , (4.12) 
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where the second equahty is derived by means of (14. 7p and (|2.44l) . (|4.9p . Thanks to p.l2p 
this means that for x — ^ oo we proved that 

^(^) = (^'S^.e^-"^^-), . . . (-))«:r' + . . . , (4.13) 

^(x) = -^v„-'[{pl:y\~^-^-\ + - . . . (4.14) 

Asymptotic values of ip{x) and ipi^) depend on the choice of matrices V and D' (cf. discus- 
sion after (|2.43p ). while asymptotics (14. lip and (|4.12p are invariant. In fact we can give a 
more detailed description of this asymptotic behavior. For this sake let us divide sectors (t„ 
in two subsectors 



o-'j = {a; : Kn+Nt-i{x) > Kn^i{x) and Kn^i{x) > Kn+N^i^)}^ 
a',[ = {a; : i^„+Ar,(a;) > Kn^i{x) and Kn{x) > Kn+N^{x)}. 



(4.15) 



It is easy to see that 

a,,^a'„Ua';Ur'„, (4.16) 
where is the ray belonging to the sector (T„ 

= {a; : Kn-i{x) = Kn+^A^) and {K„+Nt - Kn~i)x2 < 0}. (4-17) 
In analogy to (j2.32l) we can determine it by means of the directing vector 

y'n = (l^l+N, - 'il-l,Kn-l - Kn+jvj- (4.18) 

In the case Na = 1, i.e., A^f, = TV — 1 by (12.21) . these definitions are senseless as it follows 
from (|2.17p . So we will consider this case below separately. Now in analogy to Definition [2TT] 
we introduce 

Definition 4.2 We say that x ^ oo in the closed suhsector (notation x — ^ oo), or in 

the closed suhsector cr" (notation x — ^ oo), n ~ 1, . . .TV, if there exist nonnegative a and 
(3 such that a + j3 ^ +oo and 

for X oo : x — ayn-i — Py'^ is bounded, 

for X — ^ oo : x — ay'^^ — /3y„ is hounded. 

Then we have the following lemma. 

Lemma 4.1 1. If x oo then Kn-i{x) > Km{x) for any m = n + N/,, ... ,n + M —2 
and Km{x) > Kn+Nh-ii^) for any m — n, . . . ,n + Nf, — 2; 

2. if x oo then Kn+Nbi^) > Km{x) for any m = n + Nh + 1, . . . ,n + M —1 and 
Kn{x) > Km{x) for any m — n + l,...,n + Ni, — I. 



Proof. Thanks to Definition 14.21 Kn-i{x) — Km[x) = Q!.9ri-i,m,n-i + Pgn-i,m.n-i and 
Krn{x) - Kn+N,-iix) = a5m,«-i,ri-i + Pgm,n+Nt-i,n~i, whcrc gi^jn,n IS defined in (j2.19p 
and 'gi^m.n is gi,m,n with substitution Ni, ^ Ni, + 1. Then the first statement follows by 
lemma [1?TI In the same way we get K„+Nbix) - Km{x) = agn~i,m,,n-i + Pgn,m,,n and 
Km{x) — Kn{x) — agrn,n,n-i + Pgm.n,n and the sccoud statement results again from Lemma 
OB 
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Thanks to this Lemma we get from (I4.1ip 



oo, 



\ / l,m 



OO, 



(4.19) 
(4.20) 



for any n = 1, . . . , TV and m = n, . . . ,n + Nf, — 1. In the exceptional case Na = 1 it is easy to 
see that thanks to (12. 17^ both these equahties coincide, and we can use any of them for the 
asymptotic behavior in the whole sector cr„. Correspondingly, relations for the asymptotic 
values of ^k,^{x) follow from (|4.7p and (|4.12p . where it is necessary to take into account 
that images of subsectors ()4.15p are different from. o''J^_^j\j , cf. f)4.5|) . Exactly, we 

get decomposition of sectors an different from (|4.16p : 



K+atJ UK%J, (4.21) 
where square brackets denotes by (|4.4p substitution Na ^ Ni,. So that by (|4.15p 

Wn+Na] = ■ Kn+Nb~l{x) > Kn-l{x) and Kn{x) > Kn+Ni^lix)}, 



[a'^] = {x : Kn+Nb-'iix) > Kn{x) and Kn{x) > Kn+NA^)}- 



(4.22) 



The same remark as above must be given in the case A^^ = 1. 

In the following we need the asymptotics of $(a;,iK„i) and '^k,„{x) for fixed m in de- 
pendence on n, while in formulas above we have its asymptotic for fixed n in dependence 
of running m. These asymptotics can be, easily, obtained from (|4.19l) . (|4.20p and analogous 
formulas for '^i^^(x), and we can state the following theorem. 

Theorem 4.2 The leading asymptotic behavior of the discrete values of the Jost solutions 
$(a:,k) and 5'(a::,k) is given by 



^{x,iK,n) = 



V. \ / l.m 



X 
X 



OO, n € [m + 1, m + A'a], 

oo, n e [m + TVa + l,m + A/], 

OO, n e [m + TVa + 1- »7i + A/], 



(4.23) 



i-^^^ix) = 



X > OO, 



e[m + Na + l,m + J\f], 



= < 



(n) 
Pn+Ni-1 



iriA^g-if„+„,_i(.)^ 2; OO, n£[m + l,m + Na] 



(VVn 



(n) 
Pn-1 



X OO , n € [m + l,m + Jya\ 



(4.24) 
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We see that, for n = m + 1, . . . ,m + Na, the Jost solution has the same 

asymptotics along the directions on the a;-plane obtained by considering the union of the 
corresponding sectors cr„, that is, thanks to (|2.31[) . along the following directions 



u ^ 



n— m+1 



(4.25) 

i 

{K,n+Na{x) - K,n{x) > -00} n {K„i+nA^) " Kj^ix) > -Oo}, Na < Nb, 

{Km+N^{x) - Km{x) > -00}, Na = Nb, 

\{K„,{x) - K^+N^ (x) > -00} n {K„,{x) - K„,+N, (x) > -00}, iVa > ^f,, 

where, say, condition Km+N^i^) ~ Km{x) > —00 means that x — > 00 in such way that this 
difference is bounded from below. Geometrically the set in (|4.25p essentially depends on 
relation between Na and Nb- Thus we have, in the first line, a sector with a sharp angle at 
the vertex (being the intersection of two halfplanes), in the second line, a halfplane and, in 
the third line, a sector with a blunt angle at the vertex. 



5 Annihilators 

In the theory of the one-dimensional Sturm-Liouville operator it is well known that the dis- 
crete value of the Jost solution corresponding to the one soliton potential (i.e., the potential 
in (|1.4I) with ni = — K2) decays exponentially on the xi-axis as e^'^^^^l. The Jost solution of 
the heat equation, obviously, cannot obey analogous property on the whole x-plane because 
of the X2-dependence of the exponential factor (see, e.g., (|4.6p ). 

In order to find a hint for a proper formulation in the case of the heat equation, let us 
consider the trivial example of the one soliton potential, i.e., the case Na — Nb — 1 and 
arbitrary m and K2 obeying ([231). By (|2l^ and (|2lll) we have t{x) = de^i(^) + e^^(^), 
where d is a constant, that thanks to (12.71) gives (|1.4[) . Then from (|4.6|) we have 

<^{X, IK2) = -d^X, IKi) = ^^K,i.) ^ ^K,(.) ' (5.1) 

and for the asymptotics of the discrete value of the Jost solution 

\r: %? (5-2) 

>, X — > 00, 

where closed sectors <7i and tij are given by 

5T = {Ki{x) — K2{x) > -00} = {xi + {ki + K2)x2 < +00}, 
^2 = {K2ix) - Ki{x) > -00} = {xi + (ki -I- K2)a;2 > -00}, 

where we understand inequalities here like in definition (14.11) . Let, now, denote the scalar 
product of the two-dimensional vectors x = {xi,X2) and q = (91,(72) as 

qx ^ qixi + q2X2, (5.4) 

and let 

qnin = 92 - (Km + ^n)?! + KmK„ (5.5) 
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Figure 2: One-dimensional case 

for any m,n G Z. Then, thanks to ()5.2p . we get for the asymptotic behavior of the functions 
e~i^^(x,ini) and e~^'^$(x, 1^2) 



^1 eJip[-{qi ~ K2){xi + {ki + K2)X2 - qi2X2], X OO, 

1 exp[-{qi ~ Ki){xi + {ki + K2)x2 - qi2X2], X 00, 



(5.6) 



and we deduce that they are bounded for all x iff qi2 — and ki < gi < K2, this on the 
segment of the line qi2 = on the q-plane cut off by the parabola 52 = 9i (see Fig. [2]). 

This trivial result has unexpected generalization in the case Na ^ Nf,. Precisely, for 
Nb > Na > 1, the set of the g-plane, where the discrete values $(a;, inm) of the Jost solution 
multiplied by the exponent e~'^^ are bounded, becomes a polygon. As regards the dual Jost 
solution, the functions e'^^'i' k,^{x) are bounded inside a polygon for Na > > 1. 

Let, then, consider the case Ni, > Na > 1 and let us introduce the angular sectors 

Pk{q) = 0{qk,k-Na,)0{qk+N^,k) (5.7) 
where the 0's are the step functions and 

qmn — (^m ^n)qmn ^ (^■^) 

and let us call 

the points of intersection of the line qmn — with the parabola q2 = q\ in the g-plane. 

The angular sector pfc(g) has vertex on the parabola at the point Qk. The ray bordering 
the sector from the left (looking from the vertex inside the sector) is crossing the parabola 
at the point Qk-N^ and the ray bordering the sector from the right (looking from the vertex 
inside the sector) is crossing the parabola at the point Qk+Na ■ If add to the parabola 
the point at infinity and, then, we consider it as a closed curve, as k increases the sector 
move along the parabola in the anticlockwise direction. 

Then, we can prove the following Lemmas. 
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Lemma 5.1 The region of the q-plane defined by the characteristic function 



m+Nb 
k=m+Na 

is a polygon included in the parabolic region q2 > q\ with a vertex at the point Qm- More 
precisely this polygon is included in the polygon with vertices the points Qm- Qm+i, 
Qm+Nb belonging to the parabola, coincides with this polygon in the case Na — 1 and belongs 
to the strip 

min Hi < qi < max ki. (5-11) 

m<l<m+Nb m<l<m+Nb 

Proof. The angular sector pk has a.t k — m + Na the left (looking from the vertex inside 
the sector) ray crossing the parabola at the point Qm- Increasing k the sector rotates along 
the parabola in the anticlockwise direction up to the sector pm+ tv^ i that has right (looking 
from the vertex inside the sector) ray crossing the parabola just at the point Qm- Therefore 
the sectors pk for k running in the interval k ~ m + Na, . . .m + N^ cover a common region, 
which is a polygon inside the parabola with a vertex at the point Qm- For Na = 1 each 
sector Pk has right ray coinciding with the left ray of the sector pk+i and, therefore, the sides 
of the polygon are just the intersections of the rays bordering the angular sectors pk with 
the region inside the parabola and the vertices of the polygon are the points Qm, Qm+i, 
. . . , Qm+ Nb belonging to the parabola. For Na > 1 the point Qk+ Na intersection of the 
right ray of the sector pk with the parabola lies to the left of the left ray of the sector pk+i 
and, therefore, does not belong to the polygon and we deduce that the polygon has only 
one vertex belonging to the parabola and, precisely, the point Qm and that it belongs to the 
polygon with vertices the points Qm, Qm+i, - - - , Qm+Nb- In addition, we can deduce that 
the polygon defined by (|5.10p belongs to the strip in (|5.1ip .B 

Lemma 5.2 The polygon defined by the characteristic function 

m+Nb 

em{q)= n ^(5fe+JV„,fc) (5.12) 

k—m 

satisfies the equation 

(I Nb>2Na-l 

/ \ / \ J 7n-\-Na — 1 , ^ 

em{q)-Pm{q)i ]"[ eiq^+N^A N,<2Na~l- ' ^'''^^^ 

y k=m+Nb-Na + l 

and, therefore, belongs to the polygon Pm{q) considered in Lemma \5.1\ and is included in the 
parabolic region q2 > q\ of the q-plane. Moreover, the polygon defined by \5.12]) is not void. 

Proof. Thanks to the equality ]Xk=m'' ^(^k+N^.k) = lVk=m+Na^(^''-k~Nj obtained by 
shifting k -+ k — Na, we get that 

m.+Nb /m+Nb \ m+J\f 

W 0{qk+Na,k)=\ W e{qk+Na,k)\ n ^(5fe,fe-Afa), 

k—m \ k—m / k—m+Na 

as 6''s of the second factor are all present in the first one. Then, factorizing Pm{q) from 
the r.h.s. and shifting the running index k for the remaining 0's in the two products for 
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k ^ k — Na and fc — fc + Na-, respectively, we get 

n ^&k,k-NM n ^&k+N^,k). (5.14) 

k=m+Na / k=rn+Ni-Na + l 

If Nf, > 2Na — 1, the 0's in the two products in the r.h.s are already present in Pm{q) and, 
then we get the first equality in (I5.13p . If Nb < 2Na — 1, we have 

m+2Na.-l m+Nt, m+2Na-l 

n ^(%,fe-jvj= n ^(*,fe-jvj n ^(^k^k-wj (5.15) 

k=m+Na k=m+Na k=m+Nt + l 

m+Nh-1 / m+Na-1 \ m+Nb 

n 0($fe+jv„,fc)= n ^(%+iv„,fc) n ^(*+iv„,fc), (5.16) 

fc=m+JV6-JVa + l \k=rn+Ni,-Na + l / k=m+Na 

and, then, recalling the expression for Pm{q) in (|5.10p . (|5.7p and noticing that the second 
product in (j5.15l) and the first product in (j5.16p coincide, we get the second equality in 

In order to prove that the polygon defined by ()5.12p is not void it is sufficient to prove 
that the product in the second line of the r.h.s. of (|5.13p equals 1 at the point Qm- In fact we 
have from definitions (|2.19p and ()5.8p that qk+N^,k{Qm) = -~gk+N^,m,k+Na, which, thanks 
to Lemma 12. 1[ is greater or equal to zero only for m = fc, . . . , fc + Na (modjV), which is 
impossible for k ^ m + Nb ~ Na + I, . ■ . ,m + Na - I and Nb < 2Na - 1. ■ 

Lemma 5.3 The characteristic functions em{q) introduced in Lemma \5.2\ coincide with the 
characteristic function. 

m+Nt m+Nh-1 

em(g) = Yl ^(^k+N^^k) Yi ^{^l+N^,l+l), (5.17) 
k—m k—m 

Proof. It follows if we notice that the parameters kjv„+i, i^i, and kj+i for any value of I 
result to be ordered only in one of the two following ways: KjVa+i < ni < ni+i or < ni+i < 
KNa+i with the only exception of the case / = J\f (modA/"), where Kj\f+i < kn^^+i < I'^N'- 
Taking into account that q is inside the parabola q2 = qi, we get for any value of I that 
0iqi+N^.i)O{qi+N^,i+i) = 0{qi+Na..i)- ■ 

Now, we can prove the following Theorem. 

Theorem 5.4 1. In the case Na < N^, the function e~^^$(x, iKm), for any m = 1, . . . ,7V, 
is exponentially decreasing with respect to x for any q inside the polygon defined by the 
characteristic function 

m+Nb 

£m{q)= n ^&k+N^M), (5.18) 
k—m 

while for Na > Nb there exists no such domain on the q-plane. 

2. In the case Na > Nb the function e'?^^K^(a;), for any m = 1, . . . ,J\f , is exponentially 
decreasing with respect to x for any q inside the polygon defined by the characteristic 
function 

£m{q)= n ^&k,k+Nj, (5.19) 
k—m 

while for Na < Nb there exists no such domain on the q-plane. 
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3. These polygons are all included in the parabolic region (72 ^ ?! • 



Proof. In (|4.23p and (|4.25p we have seen that the geometry of the sector S„i = Un^m+i 
where $(a;, inm) has asymptotic behavior e^^-i^) essentially depends on the relation between 
Na and Ni,. Then, first, let Na < N}, and let us consider the directing vectors of the rays 

bordering the sectors Em, cr'n and cr" {n = m + Na + 1, ■ ■ ■ ,m + A/"), where the prod- 
uct e~'^^^{x,iKjn), according to (|4.23p and (|4.25p . has, respectively, asymptotic behavior 
^KUx)-qx^ gX„_i(:E)-g^^j^j ^K^+^^{x)-qx ^ p^j. ^j^^y Vm+N^, for < they are 

Vn-i and y^, and for ct" they are and (see (|2.32p and (|4.18p ). According to Definition 

I2.1l we say that x ^ 00 in the sector Em if a; — aym — Pym+N^ is bounded while + — ?• 00 
and both a and /3 are bounded from below. On the other side, it is easy to see that 

K„,{x) - <lAx=ay,^+l3y,^+r^^ " aq„^+N^,7y^ ~ Pq,n+N,,7n- (5.20) 

and we deduce that e~'^^^{x,iKm,) is bounded in the sector E„i iff qm+Nb,rn < and 
qm+Na,m > 0. By Definition 14.21 we say that x — > oo in the sectors a'^ and a'^ if x — 
ctyn-i ~ Py'n and X — ay'^ — /3y„, respectively, are bounded while + — ?► oo and both a 
and /3 are bounded from below. Since we have 

(5.21) 

and 

Kn+N,{x) - <ix\^^ay'^+l3y,^ " aqn+N^,n~l + Pqn+N^.n, (5.22) 

we conclude that e~^^$(a;, iKm) is bounded in the sectors and cr" iff q„+jVi,-i,n-i, 
Qn+Nb.n-i and qn+Nt,n are less or equal to zero. 

Thus the condition that the product e^''^^{x^iKm) decays asymptotically in any direc- 
tion on the x-plane is equivalent to condition that q belongs to the set on the q-plane given 
by characteristic function 

m+Nt m+Nb-1 

i,n{q) = W 0{qk+Na,k) W 0iqi+NaJ+i), (5.23) 

k—m k—l 

where we took into account that indices are defined by mod A/" and that qm,n — —qn.m (see 
and dS^l). 

Thanks to the four Lemmas above, the first part of the statement [T] of the Theorem for 
Na < Nb is proved. 

In the case Na = Ni, we have by (I2.32p that ym+Na — ~ym, m = 1, . . . , Na- Thus it is 
enough to consider the asymptotic behavior along the ray (see Definition [OJ. We choose 
ym as directing vector and by (|5.20p for (3 = we have to consider limits in both directions: 
a — > ±oo. It is clear that the r.h.s. of ()5.20p . independently of the sign of qm,m+Ni, cannot 
decay in both these directions and this r.h.s is bounded only if qm,m+Nb — 0, that, thanks 
to (|5.5p and (|5.8I) . gives a line and not a domain on the q-plane. Analogously, in the case 
Nb < Na the directing vectors of the sector in the third line of (|4.25p are the same like in 
(|5.20l) . but conditions on the behavior of a and /3 are different. Together with the limit 
a —¥ +00, /3 — -|-oo, we have to take into account sectors a — > -\-oo, /? — — oo and a — > ~oo, 
P — > -|-oo. Thus we see that also in this case it is impossible to give conditions on qm,m+Nb 
and qm.m+Na that always make the limit of the r.h.s. of (j5.20p to be equal to —oo. This proves 
statement [T] of the theorem. Statement [2] for e''^'i>K.^{x) follows from (|4.7p . In particular. 
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existence of the polygon only in the case Nj, < Na results from above and (|4.4I) . The third 
statement of the theorem is proved in Lemma 15.21 ■ 

Thus Theorem 15.41 is proved. Using (|2.3p . (|2.17p . and (|5.8p we can rewrite (|5.18p more 
explicitly: 



£m{q) 



J]^ 6'((7fc,fc+AfJ J]^ 6'(-gfe,A:+ArJ, 1 < m < Na, 

e— m k—1 

Y[0{qk,k+NjY[0{^qk^k+N,) H Na + 1 < m < Nt, (5.24) 

Nb + l<m<J\f. 



k—rn 
N. 



k=l 



k = Na. + l 



m-Na 



^ k—m—Nf, 



k=l 



Fig. [3] shows, for a the potential u{x) with Na = 4, 7Vb = 6, the polygon with charac- 
teristic function Smiq) on the q-plane in the case m — 5 and the union of all polygons for 
m = 1,...,10. 




Figure 3: Na — 4:, Nb — 6, JV — 10, polygon for to = 5 and union of polygons for m = 
1,...,10 



6 Concluding remarks 

Theorem 15.41 shows that the extended L-operator 



L{x, x'; q) = {-d^ - ^2 + {d-x, + qif}S{x - x') ~ u{x)S{x - x') 



(6.1) 



is not invertible in the case where u{x) is {Na, iVb)-soliton potential. If Na < Ni, it has right 
annihilators, and in the case of Nb < Na the left ones. 

In fact, in the case Na < Nb, if we introduce an operator Hm with kernel 



Hm{x, x; q) = £m(q)e '^'^^{x, Km)i){x ; q) 



(6.2) 
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where ip{x';q) is any arbitrary self-adjoint function, bounded in x' we have 

L{q)H„,{q) = 0. 



(6.3) 



Analogously, for the case Ni, < Na, we have 

K„MLiq) = 0, (6.4) 

with 

K^{x,x';q) = £™(q)(p(x; (7)e«"^'*«„(x'), (6.5) 

where ip{x] q) is any arbitrary self-adjoint function, bounded in x inside the polygon defined 
by e,n{q). 

We note that the total Green's function G(x,x',k) of the operator C can be defined 
(see [m |3]) as the value of the kernel M{x, x'\ q) at q^ = kg., 52 = k| — for a complex 
spectral parameter k — ks^-l-ikg. These values of q lie outside the parabola q2 = qf and 
therefore outside the polygon and touch it only at the vertices. The Green's function then 
exists for any k but it is singular at the points k = iKm corresponding to the vertices of the 
polygon touching the parabola. 
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Appendix 

A Regularity of potential and Jost solutions 

We defined the Jost solutions 3>(a;, k) in a symmetric form in p.3p . so that 

X(x,k) = $(x,k)e*''^i+'''^= (A.l) 

became a polynomial in the spectral parameter k with the higher power term given by 
(ik)^"-. Therefore, it can be determined by giving its values in Nt points. If these points 
are chosen to be k = iK„ with n belonging to the subset J = {ni, • • • , fiAff,} of numbers in 
the interval {1, . . . , A/"}, these values can be determined by a combined use of the analyticity 
properties of x(a;,k) and condition $(iK)2? = in (|3.9p . 
Let 

A(k)= J](ik+K„). (A.2) 

nS./ 

Then the ratio x(a;, k)A^^(k) is normalized to 1 at infinity and the poles at the points 
\i = iHn with n G J are the only departures from analyticity. We have, therefore, 

X(a^:k) ^ -, \p xix,iKn) . . „^ 

A(k) ^;^^(k-*«„) A' (*..„)■ ^ 
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This gives 

^{x,iKm) = ^{x,iKm), m£j, (A. 4) 

, m ^ J, (A. 5) 



Then condition ^(in) 2? = is equivalent to 
Comparing with (|4.6p we get for the r-function (up to the standard similarity property) 



= - ^ A(iK„)e'^" ■ (A.6) 



r{x) ^ det ( -^ ^ , ) . (A.7) 

^ l=l....N^ 

By inserting this r-function into p.7p we get the soliton solution and different choices of the 
set J correspond to different equivalent formulation of the soliton solution, in accordance 
with the result presented in [TB]. 

We want to get a condition of regularity for the potential equivalent to (|2.18p . but 
involving only the Na x A';, submatrix c?(„) introduced in p.44[) and not the full J\f xNi, 
matrix T). 

This, of course, can be obtained by pure algebraic methods. However, we want to get 
this result by using the freedom in expressing the t(x) function offered by (jA.7p . showing 
in this way its usefulness. 

To the end of getting this result we consider the evolution with respect to a multivalue 
time describing the evolution of the entire KPII hierarchy. Then, the exponents in the 
different formulations of (|A.7|) are independent and the requirement that their coefficients 
are grater or equal to zero in one formulation implies the same requirement in all other 
formulations. 

We can, then, state the following theorem. 

Theorem A.l Let T) he a J\f xNb matrix of the form 

V = ttB, with 6 = ^ , (A.8) 

where d is a real Na x Ah matrix and E(, is the unit N(, x A^ matrix and where the J\f x Af 
matrix 

^ — ll^mn ||m,n=l 7 ""mn — ^■!T{m),n = ^m,-!T~^ (n) i i-^-^) 

performs a permutation of the M rows of the matrix on the right according to the transfor- 
mation 

tt: (l,...,A/-)^(7r(l),...,^(AA). (A.IO) 

Then the matrix V has all its maximal minors non negative if and only if the Na x Nf, matrix 

dfcr = fc = l,...,7V,, r = l,...,Ar,, (A.ll) 

where {ji^ j2, ■ ■ ■ ^ Jn^} ^■s a permutation of {1,2, Na} such that 

< < ■ ■ ■< TT-^O'jvJ, (A.12) 
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and {h, Z2, . . . , In,,} *s a permutation of {1, 2, . . . , N^} such that 

T^-\N, + h) < 7T-\Na + /2) < ■ • • < 7r-\Na + InJ, (A.13) 

is, according to the definition in U9f . totally non negative, i.e. has all its minors non 
negative. 

Proof. From (|A.8|) and (|A.9|) we have 

(A.14) 

Let, now, J be the set 

J^{7r-\Na + l),...,7r-\N-)} (A.15) 

and let 

CJ={7:-\l),...,7:-\Na) (A.16) 

be its complement in the set {1, . . . ,Af}. 
Then, from (IA.7I) we have 



(*«:7r-i(Ar„+o) 

1=1 



Na 



= - ^A(iK,-i(,))e^~'-^o)(i,-/-. (A.17) 
For the r function we obtain, therefore, 

r(x)=detUM'+*E7 wTT V^^-^H " (^'l^) 

\ / l,l' = l,...,Nt 

Let us re-order the rows in the matrix in the determinant from Z = 1, 2, . . . , A'^t, to a permu- 
tation li,l2, . . . ,lNb such that 

KT,~i{N^+h) < Kt,-1(n^+12) < ■ ■ ■ < '«7r-i(JV„+ijvJ (A. 19) 

and let us re-order the sum over j from j — 1,2,..., Na to a permutation ji , j2, . . . , jAf„ such 
that 

^^-Hn) < '*T-i02) < • • • < i^-n-HjNj- (A. 20) 

Then, up to an unessential sign, the t function can be rewritten as 

t{x) = dot ^Er,l' + 

+ 1 ( a^!;::;i:;;,) ) „; •^■^^ 

where 

^^r./' = Si^,i> (A.22) 



23 



A,fe = . (A.23) 

If we multiply from the right the matrix in the determinant by we get up to a sign 

t{x) = det(£;A,, + T) (A.24) 

where {r,s = 1, . . . , Nb) 



(A.25) 



By using repeatedly the Binet-Cauchy formula, we have 
and, then. 



^ :!: :'! = E n 



X (n AK-.u.„,). ''■-■■..-) A ( )«( ) " 

Now, if we impose that the coefficients of the exponents are non negative, we get 

y Ki, K2, ■ ■ ■ , Kn J Y ri, r2, . . . , r„ y 

Since 

ri,r2,...,r„ \ i<i<m<n 



A 



n 



z,m— 1 

(A.29) 

the denominator in (jA.29|) cancel with 11™=! ^(*'*7r-i(jfc„)) in (|A.28|) . and, after reversing 



the order of the r^'s, we have 



nA(-.-M..)))Ar-g;;;;;;; J>o. (a.3o) 



Finally, since 



zA'(iK^-l(Ar^+;^^)) = (-1)'"" ^ |-'«7r-i(A'a+ir„J + '^^-MJVa+irjl > (A. 31) 
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we get the condition 

f ^n, K-u. . . , fci \ g_ 

\ T'\_,T2, ■ ■ ■ iTn J 

We conclude that the matrix 

dkr = {-lf+''-d^(o,,^^,_,,),W, k = l,...,N,, r = l,...,Nt, (A.33) 

is totally non negative. ■ 

Prom this theorem the following Corollary follows easily. 

Corollary A.l A Na x Nb matrix d is totally non negative, i.e. has all its minors non 
negative, if and only if the J\f xNb matrix 

where _ 

dji = {-iy+'dN^+i.j,i, j = l,...,Na, l = l,...,Nb, (A.35) 

and Eb is the unit Nb x Nb matrix, has all its maximal minors non negative.M 
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